CONVERGENCE OF AN ADAPTIVE APPROXIMATION 
SCHEME FOR THE WIENER PROCESS 

MATS BRODEN AND MAGNUS WIKTORSSON 

Abstract. The problem of approximating/tracking the value of a Wiener 
process is considered. The discretization points are placed at times when the 
value of the process differs from the approximation by some amount, here 
denoted by r\. It is found that the limiting difference, as 77 goes to 0, between 
the approximation and the value of the process normalized with rj converges 
in distribution to a triangularly distributed random variable. 

1. Introduction and preliminaries 

An adaptive approximation scheme of the Wiener process is considered. The 
Oh ■ discretization points are placed at times when the value of the true process differs 

(-h from the approximation by some amount, here denoted by r\. This can be seen as a 

control problem where we want to track the true value of the process with our ap- 
proximation, and where both the process and its approximation are fully observable. 
The approximation strategy presented here may be feasible when discretization is 
associated with some cost that should be kept low. Examples of related problems 
is that of discrete time hedging of derivative contracts in financial markets (see e.g. 
P*" ' Geiss and Geissl (|20f)fir i1 and c ertain space-time discret i zation schemes of stochastic 



differential equations (see e.g. Milstein and Tretvakov (199^)) 



OO 

\ Let X be a diffusion process defined by X t = aW t , where W denotes a one di- 

■ mensional standard Wiener process. Define, for some 77 > 0, a sequence of stopping 
^ ! times {?-}i>o by 

O: % +1 =M{t>%\ \X t -X t n\=T,}, 

where = 0. The components of the sequence t v may be seen as epochs of the 
renewal process N n defined by = sup{i : < t}. Furthermore, let the sequence 

■ { T i}i>i °f interarrival times be defined by t£ = — and define the renewal- 
reward process (p by ip^ := X^i 7 ^- The process X n may also be seen as a 
renewal-reward process, but with a reward that takes the values —77 and 77 with 
equal probability. 

The aim of this work is to investigate the asymptotic behavior of (X t — X^ v ) / 7/ 
as 77 approaches 0. It will be seen that this quantity converges, pointwise for each 
t > 0, in distribution to a stochastic variable which is triagularly distributed. 

Before we end this section we will state some resluts regarding barrier crossings 
and renewal processes. The main result is presented in Section [21 In Section [3] we 
perform a simulation study and investigate the transition to the limiting distribu- 
tion. 
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1.1. The Wiener process with two absorbing barriers. Since the components 
of the sequence {T^}i>i are independent and identically distributed, we will let t 71 
denote a stochastic variable with the same properties as these t^'s, and which may 
be characterized by 

r^mi{t>0\\X t \=r ] }. 

Now. consider the process X absorbed in —77 and ry, that is X tA ^. The transition 
densit y of this process, from Xq — 0, may be represented by (see ICox and Miller 
dl965h ) 



OO 1 

(1) P "(t,x) = ^4 e -K^f) 2 * S i 



kir\ . ( k-K{x + 77) 



for all (t, x) G (0, 00) x [—77, rj\. This tr ansition density may al so be expressed as an 
infinite sum over Gaussian kernels (see Cox and Miller ( 19651 )) 



r-^ 1 / (x-4fc„) 2 (x-2^+4fc^) 2 

(2) p"(t,x)= -1^7 (e-^^-e 

k=-oo V2tt<jH V 

for all (t,x) G [0,oo) x [—77,77]. 

Lemma 1. The integral o/p 1 (7j,a;) 

a) with respect to t over the interval [a, b] C [0,oo) may be represented as 

for all x G [—7/, 77] . 

b) with respect to x over the interval [a, b] C [—77,77] may be represented as 

(3) f a P<\t,x)dx= ± / a 6 ^- |( ^ )2t ^(^)sin(^|^)^ 

/or /; G (0, 00), or as 

C f b 1 / (x-4fci7) 2 (x-2 + 4fcr,) 2 \ 

(4) / p q {t,x)dx= \ / e i^^-e ^ Wfe, 
J a k=-oo^ a v2ira 2 t \ / 

for all t G [0, 00). 

Proof, a) Define the functions gj^ and G^by 

g-K^f) * /fc7T\ . /fc7T(.T + 77) 

and G£(i, ar) = £™ =1 9k (*> x ) thcn lim «t°° G£(t, z) = p"(i, a:). Since 5 f (t, 0) > 
it follows that 

0<G^,0)<G^ +1 (t,0), 
and consequently by Lebesgues monotone convergence thorem 



/ lim G n (t,0)dt= lim / G n (t,0)di. 

Jo "too «tooJ a 



Extending the integral we get 



lim / <?„(*, 0)di < lim / G n (t,0)dt. 

"too J a ntoo ,/ 

Moving the integral inside of the sum in G„(i,0) and performing the integra- 
tion over R + we get the sum lim^oo J2k=i 8/(fc 2 7r 2 er 2 ) = 4/(3cr 2 ), and hence 
linintoo J G n (t,0)dt < 00. Since | sin(fc7r/2) sin(fc7r(a; + rj)/(2rj)) \ < 1 it holds 
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that \g%(t,x)\ < g£ (t,0) which implies that |G£(t,x)| < 0). Since G£(t,0) is 

bounded by limn-j-oo G^(t,0) the function G%(t,x) is dominated by the integrablc 
function lim^oo G^(t, 0) and by the dominated convergence theorem it follows that 



r-b r-b 

I lim G n (t, x)dt = lim / G n (t,x)dt. 

J a "too "tooj a 



Moving the integral inside of the sum on the right hand side the claim is proved. 

b) Eqn. From the proof of a) we know that G^(t,x) < p v (t,0) = lim n foo 
G n (t, 0) which is bounded for every t > 0. Since the set [a, b] is bounded (i.e. 
[a,b] C [ — 77, 77]), the claim now follows from the bounded convergence thorcm. 

Eqn. (J?]) Define the functions and G„ by 

(x — fc) j-j, 

fffc (*,*) = ^7== and G£(i,x) = 2 (gf kv (t,x) -gtiknit,*)), 

V2tt<jH k= _ n 
then lim^oo G^(t,x) = p v (t,x). The function G^ may be decomposed as 
G^(t,x) = G^ 1 (t,x) + G^ 2 (t,x), 



where 



and 



fc=0 



G£' 2 (*.z) = £ (fff fcl ,(t,x) -sf^+aM)- 
fc=-i 

Since each term in G^' 1 is positive and each term in G„' 2 is negative it holds that 

< G^ l (t,x) < G^+tfax) and > G%' 2 (t, x) > G%(t, x). 
The claim now follows by Lebesgues monotone convergence theorem. □ 
Lemma 2. It holds that 

1*00 

a 2 / p 1 (t,x)dt = {l-\x\) + . 
Jo 

Proof. From Lemma [T] a) we have that 

' kir(x + 1) 



p 1 (t,x)dt 



8^1. fkn\ 1 . 

k=l v 7 



The idea is to find a function that can be expressed as a series which corresponds 
to the above sum. Let Si = 1/2 and S2 = (x + l)/2, then 

2 ^2 poo 00 

/ p 1 (t,x)dt = N — sin (fc7rsi) — sin (fc7TS2) ■ 

8 7 £j* fc 

Define the function h s by 

\ 1 , s < \x\ < 1. 
The Fourier Cosine coefficients of /i s are given by 

Co = / /i s (x)dx = 1 — s , 
Jo 

f /1 1, , 1, 2 sin(7rfcs) 

afc = 2 / cos{kirx)h s (x)dx = - . 

Jo Kit 
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Applying Parseval's formula yields 



f 1 , , , , , . , sin(7rfcsi ) sin(7rfcs 2 ) 

/ h Sl (x)h S2 (x)dx = 2 V K 11 \ 2 > + (1 - Sl )(l - fla ) . 
Jo ,_ 1 7r ' s 7r ' c 



Assume that x £ [0,1], then 0<si<s 2 <l and 
2 

7T 2 /C 



oo 2 

^ — jrj sin(7rfcsi) sin(7rfcs 2 ) = 1 - ,s 2 - (1 - Si)(l - s 2 ) 



fc=i 



si(l - s 2 ) 



Thus 



fc-7r\ / kiv(x + 1) 
sin 1 



/ V 2 

= (l-x). 

Repeating the argument with x £ [—1,0] yields the result. □ 

One important property in the theory of renewal processes is that of direct 
Riemann integrability of a function. A function function H(-) is said to be directly 
Riemann integrable over [0,oo) if for any h > 0, the normalized sums 

oo oo 

h > inf H(nh — 5) and h > sup H(nh — 5), 

^0<S<h ^0<8<h 
n—1 n—1 — — 



conve rge to a common finite limit as h 4- (see chapter 4.4 in Dalev and Vere-JonesI 

diaii)). 

Lemma 3. TTie function p 1 ^, a;) is directly Riemann integrable with respect to t 
for each x £ [—1,1]. 

Proof. We will start by considering the case when x = 0. The function p 1 ^, 0) 
is directly Riemann integrable if p 1 (t, 0) is nonegative, monotonically de creasing 
and Lebesgue integrable (see chapter 4.4 in Dalev and Vere-JonesI ( 19881 )). Since 



each term in the representation ((3]) is nonegative and monotonically decreasing for 
x = so is p 1 (t, 0), and by Lemma [2] the integral of p 1 (t, 0) over [0, oo] is given by 
/ °°p 1 (t, 0)dt = 1 and thus p x {t, 0) is Lebesgue integrable which proves that p 1 (t, 0) 
is directly Riemann integrable. 

Next let x £ [—1, 1] \ {0}. The function p(t, x) is directly Riemann integrable 
with respect to t if p 1 (t, x) > 0, p(t,x) is uniformly continuous in t and bounded 
from above by a monotonic ally decreasing integrable function (sec chapter 4.4 in 
Dalev and Vere- Jon es (1988)). Since p(t, x) is a probability distribution for each t it 



is clear that p(t, x) > 0. To show uniform continuity we will split the interval [0, oo) 
into two parts, say [0, 1] and [l,oo), and show that p x (t,x) is uniformly continuous 
on each part. For the interval [0, 1] we will use the representation Let g^ and 
be defined as in the proof of Lemma [TJ It is clear that each gjf is uniformly 
continuous in t and thus also G„ is uniformly continuous for each n < oo. If we 
can shown that G„ (i, x) for each x £ [— 1, 1] \ {0} converges uniformly with respect 
to t over [0, 1] as n f oo, then also the limit p(t,x) will be uniformly continuous. 
Rewrite G% as G%(t, x) = ££ =0 g%(t, x) where g§(t, x) = gfi(t, x) - g§{t, x) and 

9k x ) = 94k(t, x ) ~ 92-4k(t> x ) + 9 G 4fc0> x ) ~ 92+4k{t, x )i for fc > 1. 
According to Wcicrstrass M-test, if there is a series of constants Mh such that 
J2kLo^ k * s convergent and \g%(t,x)\ < for all t £ [0,1] then G^ converges 
uniformly in [0,1] as n f oo. The functions gk{t,x) attains its maximum at t = 
(x - k) 2 /a 2 A 1 for t £ [0,1], and thus g k (t,x) < g k ((x - k) 2 /a 2 A l,x). The 
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function go(x 2 /<r 2 A l,x) is bounded and it is easily seen that the functions g^ 
may be bounded by C/(l + k 2 ), for some bounded constant C, and which is clearly 
convergent. Hence, for each x G [—1,1] \ {0}, p{-,x) is uniformly continuous in 
[0, 1]. To show uniform continuity in [l,oo) we will use the representation ([3]). Let 
t > 1, then 



-t\- k "^ S 



e 8 |e 
fe=i 

V 82 k 2 a 2 v 2 J_ 
fc=i 

where we used the inequalites e -2 ' < y~ 2 and le -2 ' — 1| < y which holds for y > 0. 
Hence for every e > we may chose S such that 5 < 4<7 2 e/3 which holds for every t 
in [1, oo). Hence p l {-, x) is also uniformly continuous in [1, oo), which together with 
the previous result yields that p(-, x) is uniformly continuous in [0, oo). In the proof 
of Lemma Q] we showed that p(t,x) < p(t,0), and that p(t,0) is a monotonically 
decreasing Lebesgue integrable function. Hence, p(t, x) is also directly Riemann 
intcgrable with respect to t for x G [— 1, 1] \ {0}, which together with the first result 
of this proof yeilds that p(t,x) is directly Riemann integrable for x G [—1, 1]. □ 

The next two lemmas regards properties of the random variable t v defined earlier 
in this section. Let F T r, denote the distribution function of t v . Lemma [5] states 
that that r 17 has a density, which we will denote by f T r, . 

Lemma 4. The expectation of is given by E[r n ] = ij 2 /a 2 . 

Proof. Let g(x ) = E[t v ], where xq denotes the initial point of t he process. The 
functi on g satisfies the following ordinary differential equation (see ICox and Miller 
(11965^ 



2 ^2 

y^( a; o) = -l, mi (-77) = m x {rf) = 0. 

The solution to this problem, with xq = 0, is given by g(0) = rj 2 /a 2 , as was to be 
shown. □ 

Lemma 5. The random variable t 11 has a density, denoted by f T v, that may be 
represented as 



f Tn (t) = / I (77 + Akm)e "2^' - (jj + 2 - Akrf)e 



I (T7 + 4fcTj) 2 _ (7j + 2-4fcTj) 2 
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(77 — 4krj)e 2^* — (77 — 2 + 4krj)e 



for all t G [0, 00). 

Proof. In this proof we will use the representation Let gj? and G„ be defined 
as in the proof of Lemma [1] By the use of Lemma [1] for t G [0, 00) 

P(r*<t) = l- / (gf kv (t,x)-gt4 kv (t,x))dx. 

k=-oo J - , l 

If each term in the sum above is diffcrcntiablc on [0, 00) and 
(5) X! a / (9?k v (^ x ) ~ 92-4k v {t,x))dx 
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converges uniformly on [0, oo) then 

j t P(T r ><t) = - J2 j t J (afk^^-gtik^^))^- 



fc=-c 



Calculating the integral and differentiating with respect to t we get for each term 
in © 

d f n 

(sSc,<X x ) " 92-4k V (^ x )) dx 



dt 



1 / , , . (v + ikr,) 2 ( n + 2-4k n ) 2 



(77 — Akrj)e aS^t — (77 — 2 + Akrj)e 2^ 



The maximum of the function e 2 <*' J t /t 3 / 2 in [0, 00) is attained at t = (a 
fc) 2 /(3cr 2 ). For the first term in the expression above we get that 

(ri + &krj) (v+iky) 2 f 3\ 3 ^ 2 a 2 e~i 1 
,; -_ e i^T- < 1 > 



2^/ 2 V2^2 "V 2 / V 2 V^ (1+4A:) 2 ' 

which may be bounded by C/(l + fc 2 ), where C is a bounded constant. In a similar 
manner it can be shown that the rest of the terms in Q may also be bounded by 
C/(l + k 2 ), and thus 

AC 



(7) 



d /"' 



< 



1 + k 2 



Since J2T=-oo 4C/ (1 + ^ 2 ) is a convergent series by Wierstrass M-test the sum 
converges uniformly on [0,oo), and hence, the density, f T v, may be represented by 
the sum ©. Since the terms in the sum of (JSJ) could be bounded by 4C/(1 + k 2 ) 
we have that \f rn (t)\ < 4C^^L_ 00 1/(1 + k 2 ) < 00 which shows that f T n(t) is 
bounded in [0,oo). □ 

1.2. Renewal processes. In this paragraph we will focus on a renewal process de- 
noted by N with idenpendent and identically distributed interarrival times {ri}i>i. 
Define the renewal function M by M t = E[N t ], and let /i denote the mean time 
between renewals, that is fi — E[ri], which holds for alH > 1. Next, we will state 
the key renewal theorem that will be needed later on. 

Lemma 6 (Key renewal theorem) . IfH(-) is a directly Riemann-integr able function 
then 

lim / Hit - x)dM(x) = - f H{x)dx . 

Proof. See e.g. Daley and Vere- Jones (|l988 ). □ 

Let F T denote the common distribution function of the stochastic variables Tj. 
Since the components of {Ti}i>\ are idependent and identically distributed the dis- 
tribution function of the sum £^ i=1 Tj may be represented by the fc-fold convolution 
of F T (here denoted F* k ), i.e. 

/ k \ 

F; k (t). 



itr') 
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Lemma 7 (Theorem 5.4 in lHevman and Sobell ( 1982 )). There exists a one-to-one 
correspondence between F T and M , and M has the representation 



M t = J2F; k (t). 



fc=l 

Under the assumption that F T has a density (here denoted f T ) we have that 

f ;\t) = j t F; k {t), 

where f* k is the fc-th convolution of the density function f T . We may now dchnc 
the renewal density m by 

, oc 

(8) m t :=-M t = Y,f: k (t)- 

k=i 

2. Main result 

In this section we state and prove the main result of this paper. To ease the 
notation in the proof we will let Z% = X r t ' — . 

Theorem 1. Fix a point t > 0, then 

where A is a stochastic variable with density function given by 

f A (z) = (l-\z\)+. 

Proof. Denote by Y t n (u) the quantity 

Y?(u)=X t -X v7 \{t-tf=u}. 

Because of the time homogeneity of the process X the following equality in distri- 
bution holds 

X t - X v7 \{t - =u} = X U \{\X S \ < 77, < s < u} . 
Consequently the density function of Y^(u) can be expressed as 

fY?(u)(y)- p{TV>u y 

The distribution function of Z% is given by 

Jo 

where 

dF t _ v n{u) = |<y(«- t)P(r" > *) + E " V? < «. N? = k)^du. 

The probability in the last term of the above expression can be rewritten as 

P(t < u, JV? = k) = P ( t- j^rV< u, j^rV < t < j^rV + r£ +1 

V 3=1 3=1 3=1 

= pt-^;<M<t-^<4i 

V 3=1 3=1 

/ / r t(«)/ T ,(z)dzd«, 
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where f T v (which exists due to Lemma [5]) is the density function of t v , and /** 
denotes the fc-th convolution of f T r,. Differentiating the above expression with 
respect to u yields 

0^ U J /#(«)/■>-> {z)dz dvj = J /;*(t - u)f T v(z)dz 

= f$(t-u)P(T*>> U ). 

This gives us that 

dF t _ v v(u) = j*(u - t)P{^ > t) + f] f*Ut - ^)P(r n > u)| du . 

Using the scaling property of the Brownian motion the following two relations are 
easily deduced 

P(r" > t) = P{t x > t/if) and -Yl<(u) = Y t l {u/if) . 

The first of the two relations above yields 

fMt) = "^(r" > t) = -^Pir 1 > t/ V 2 ) = ^fAt/v 2 ) , 

and consequently 

dF t _ v n(u) = L{u-t)P^ > t/v 2 )+jt^f:? (^) Pir 1 > V? 2 )| du. 

The relation Y^(u)jr\ = Y^(u/rf) yields 

fY^{ u )/r,{y)dF t _ v v{u) = I f Yt l (u/^)(y)dF t -^(u) , 

JU JO 

and thus 

Jo P( T > u hi ) V v J 

= P \th\ z) + j\\u/ n \z) jr ^/;f (^) du . 



fe=l 

2 i 



Now, by a change of variables (w = (t — u)/rj 1 ) 

•A) / fe=1 

Since 

877 2 4?7 2 



|pW,*)l<E 
fe=i 



k 2 a 2 n 2 3<t 2 



we have that lim ?? _j.o P v {u-, t/rj 2 ) = 0. For the second term we have using (|5J| , Lemma 
U and Lemma [5] (which applicable since x) is a directly Riemann intcgrablc 
function due to Lemma [3]) 

rt/n 2 ft \ °° f 00 
lira / p 1 (y, -j - v V #£ («) dv = a 2 / p 1 (y,w)d M . 
Jo \ T J Jo 
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Now by Lemma [2] 

lim/ z?/r; (z) = (l-|z|)+, 

as was to be shown. □ 

Remark 1. Note that the limiting distribution does not depend on a. This is 
unlike the case when discretization takes place on an equidistant grid, where a 
affects the variance of the limiting distribution. Instead, in the case of adaptive 
approximation, a is related to the expected number of discretization points. 

Remark 2. In the proof above all intcrarrival times up to the time t is used 
in order to characterize the distribution of t — tp]! . However, we belive that t — 
may be characterized by the dynamics of the process X in a small region around 
X t . which would imply that the result of Theorem Q] is a local result. From this 
and the fact that the diffusion coefficient a is scaled away in the limiting expression 
of the distribution we conjecture that Theorem [1] would hold for a larger class of 
stochastic processes such as SDE's. We plan to address this in future research. 



3. Numerical results 

In this section the transition of j z n j v as r\ goes from some large value towards 
zero is investigated. We will argue that for large values of rj the stochastic variable 
Zt jr\ is approximately normally distributed, and thus as r\ approaches zero we will 
see that f z v j v goes from the density of a normally distributed random variable to 
the density of a triangularly distributed random variable. 

A total of 50000 trajectories of the process X was simulated, over a period from 
t = to t = 0.5, with a = 1, on a time grid with 200001 equally spaced points. 
Trajectories of the approximation X^v were calculated for a number of different 
values of rj in the range [0.5, 4.0]. 

Recall, from the proof of Theorem [TJ the expression of the density 

,i/,, 2 / , x oo 

(9) fz*/ v (z) = p^t/rf, z)+ P 1 [--^AzZ f$ ^ dv ■ 

It is clear that for large values of rj it is the first term in (j9]) that is the dominant 
one. Thus, in this case the density is approximately the same as the absorbed 
Wiener process. Furthermore, since rj was assumed to be large the density of the 
absorbed Wiener process is approximately the same as the Wiener process without 
absorbing barriers. Hence, for large rj we have that 

do) h!ln(z)a ^( z jL-y 

where <f> denotes the standard normal density function. 

In Figure [T] the density of fz?/r]> at t = 0.5, as we let rj go from 4.0 to 0.5 is 
depicted. It is seen that when rj = 4.0 the distribution is quite close to the normal 
distribution. For 77 = 0.5 the distribution on the other hand is quite close to the 
triangular distribution. 

To further illustrate the transition from the normal distribution to the triangular 
distribution we measured the distance in therms of the Wasscrstein metric between 
the, from the Monte Carlo simulation, estimated distribution and these two dis- 
tributions. The distance between two distributions, with distribution functions F 
and G, in terms of the Wasserstein metric is defined by 

d w (F,G)= { \F(x) - G(x)\dx . 
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Figure 1. Left (large values of 77): kernel estimates of fz?/r]( z ) 
where 77 = 4.0 (dotted line), 77 = 3.25 (dash-dotted line) and 77 = 
2.5 (dashed line), and the Gaussian distribution (solid line). Right 
(small values of 77): kernel estimates of fz?/r)( z ) where 77 — 2.5 
(dashed line), 77 = 2.0 (dash-dotted line) and 77 = 0.5 (dotted line), 
and the triangular distribution (solid line). 



0.5 
0.4 
0.3 
0.2 
0.1 



12 3 4 

V 

Figure 2. Distance in terms of the Wasserstein metric between 
the triangular distribution and the empirical distribution (squares) , 
and the normal distribution (|10[) and the empirical distribution 
(circles). 

In Figure [5] the Wasserstein distance between the empirical distribution and the 
triangular distribution as well as the distance between the empirical distribution 
and the normal distribution (fT0|) . at t = 0.5, as a function of 77 is depicted. Note 
that in the case of the normal distribution ([TP]) not only the empirical distribution 
but also the normal distribution that we compare with is dependent of 77. It is 
seen that for 77 smaller than 1.25 the empirical distribution is relatively close to 
the triangular distribution whereas for values over 2.25 it is close to the normal 
distribution ([TO]) . For 77 in the interval (1.25,2.25) the distribution is probably 
better explained by a mixture of the two distributions. The small offset from zero 
for small values of the distance is due to the variance of the monte carlo simulation. 

From it is clear that it is possible to fix 77 and instead of letting 77 approach 
zero let t approach infinity. To capture this we have plotted the variance of jr\ as 
a function of t for a couple of different values of 77 (see Figure[3]). The constant 1/6, 
that is the value of the variance of the triangularly distributed random variable, 
is also plotted in the figure. As expected it is seen that for low values of 77 the 
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0.1 0.2 0.3 0.4 

Time 



Figure 3. The variance of Z% Ji) as a function of time where rj = 
0.50 (dotted line), 77 = 0.75 (thin dash-dotted line), r, = 1.00 (thin 
dashed line), 77 = 1.50 (thick dash-dotted line) and rj = 2.25 (thick 
dashed line), together with the function (t/0.5 2 )A(l/6) (solid line). 

limiting variance of 1/6 is attained much faster than for higher values of 77. From 
the argumentation above regarding high values of r\ it is also clear that for low 
values of t the distribution is approximately normal. Hence, the slope of the lines 
near zero is given by 1/r/ 2 , as is seen in the figure. 
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